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Abstract 

Simultaneous use of EEG and fMRI has become a powerful 
and attractive multimodal brain imaging technique in recent 
years, and has been broadly employed in neuroscience 
research as well as in clinical practice. The possibility to 
integrate the high temporal resolution of EEG with the high 
spatial resolution of fMRI constitutes one of the relevant 
advantages of using combined EEG-fMRI. Meanwhile, a 
number of challenges have to be overcome in order to 
consolidate such a technique as an independent and effective 
method to brain imaging. In particular, the artefacts which 
arise in the EEG signal, induced by varying gradient 
magnetic fields within the fMRI magnetic scanner. This 
work presents a novel methodology for waveform artefact 
modelling in order to subtract and clean up the EEG 
recordings. Implementation of our method results from the 
combination of two techniques: a non-linear low-pass filter 
approach based upon the slope adaption between 
consecutive samples of the signal (SSD); and the modelling 
proposal for the underlying gradient artefact components 
which is implemented using the sum of a set of sinusoid 
waveforms. The resulting EEG restoration obtained by our 
methodology shows to be similar to the artefact correction 
achieved by the established average artefact subtraction 
(AAS) method. Moreover, the proposed mathematical model 
for the artefact components shows to predict the variability 
of the artefact waveform over the time. 
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Introduction 

Acquisition of the electroencephalogram during 
functional magnetic resonance imaging (combined 
EEG-fMRI) is a multimodal technique for brain 
activity mapping which has achieved a wide usage in 


research and clinical purposes due to its powerful 
capability to provide new insights into the brain 
function (Villringer et al., 2010). The advent of the 
functional magnetic resonance imaging (fMRI) around 
two decades ago (Belliveau et al., 1991) allowed the 
combination of this method with 
electroencephalography (EEG) in studies of epilepsy 
(Warach et al., 1996; Seeck et al., 1998). Combined 
EEG-fMRI revealed itself not to be only an additional 
monitoring tool, but a promising technique for brain 
activity mapping. Nowadays, it has been extended to 
other types of neuroscientific studies and attracted the 
interest of several researchers and clinicians (Ritter 
and Villringer, 2006). 

According to Villringer et al. (2010), as being 
complementary methods, only simultaneous EEG- 
fMRI offers the opportunity to relate brain imaging 
modalities to actual brain events, because the high 
temporal and spatial measurements which are 
achieved by the combination of both techniques. This 
property may help to understand numerous questions 
in basic and cognitive neuroscience. In this way, 
combined EEG-fMRI have been used to investigate 
cognitive measures of memory performance, learning 
capabilities, location of epileptic foci, and 
characterization of the relationship between epileptic 
electric activity and hemodynamic response (Villringer 
et al., 2010; Ritter and Villringer, 2006). 

Parallel to the breakthroughs achieved by using 
simultaneous EEG-fMRI, technical problems 
associated with the induction of spurious voltages in 
the EEG in the magnetic fields of the fMRI scanner are 
urgent to be solved in order to consolidate and 
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broaden its applications range. That is the case of (a) 
artefacts caused by the subject movement within the 
scanner; (b) the "ballistogram" artefact, induced in the 
EEG electrodes by the pulsatile movement of the 
blood in scalp arteries within the static magnetic field 
of the fMRI equipment (Bo) (Allen et al., 1998); and (c) 
the "gradient" or "imaging acquisition" artefact, 
caused by the variation of the gradient magnetic fields 
within the fMRI scanner (Allen et al., 2000; Mulert and 
Hegerl, 2009; Ritter et al., 2010). 

According to Ritter et al. (2010), gradient artefacts 
occur in the EEG signal due to the voltage induced by 
the application of rapidly varying gradient magnetic 
fields for spatial encoding of the magnetic resonance 
signal (MR), and radiofrequency pulses (RF) for spin 
excitation in the circuit formed by the electrodes, leads, 
patient and amplifier. The waveform of the gradient 
artefact caused by one gradient pulse is approximately 
the differential waveform of the corresponding 
gradient pulse (Anami et al., 2003; Ritter et al., 2010). 
Imaging acquisition artefacts have amplitudes that can 
be several orders of magnitude higher than that of the 
neuronal EEG signal. Artefact amplitudes associated 
with the gradient switching (10 3 to 10 4 pV) are 
generally much larger than those arising from RF 
pulses (up to 10 2 pV) (Ritter et al., 2010). As discussed 
by these authors, gradient artefacts have a strong 
deterministic component because of the pre- 
programmed nature of RF and gradient switching 
sequence (EPI sequence). In addition to exceeding the 
bandwidth of the standard clinical EEG, the gradient 
artefact frequencies overlap the EEG frequency band 
in discrete harmonic frequency intervals or "frequency 
bins", whose fundamental corresponds to the inverse 
of the EPI slice time parameter (Niazy et al., 2005). 

In the literature, various techniques have been 
suggested in order to minimize the effects of the 
gradient artefacts in the EEG signal. For example, it is 
possible to reduce their magnitude at the source by 
laying out and immobilizing the EEG leads, twisting 
leads, using a bipolar electrode configuration, and 
using a head vacuum cushion. Furthermore, 
depending on the application, a periodic interleaved 
approach, whereby the MR signal acquisition is 
suspended at regular intervals, could be used as well 
(Ritter et al., 2010). Mullinger et al. (2011) also 
proposed a procedure for reduction of the magnitude 
of gradient artefacts by adjusting the subject position 
within the fMRI scanner. Regarding continuous MR 
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acquisition, dedicated software solutions have to be 
developed for gradient artefacts removal in the 
electroencephalogram (Con calves et al., 2007). Thus, 
post-processing signal methods like frequency-domain 
filtering (Hoffmann et al., 2000; Sijbers et al., 1999; 
Sijbers et al., 2000), average artefact subtraction (AAS) 
(Allen et al. 2000), principal component analysis 
(Negishi et al., 2004; Niazy et al., 2005), independent 
component analysis (Mantini et al., 2007), and spatial 
filtering (Brookes et al., 2007) have been proposed and 
employed for artefact correction. 

Amongst these techniques, the average artefact 
subtraction (AAS) methodology, the most popular for 
imaging artefact correction, consists of the calculation 
of an average artefact waveform or template over a 
fixed number of samples, which is then subtracted 
from the EEG signal at each sample (Allen et al., 2000). 
However, although the AAS allows a substantial 
removal of the gradient artefact, it is assumed that the 
arterfact template has a constant morphology over the 
time and, therefore, the variability of the artefact 
waveform is excluded from consideration. Thereby, 
after the subtraction is performed, artefact residuals 
arise in the signal and further processing methods 
must be carried out in order to achieve a satisfactory 
signal correction. Moreover, effectiveness of the AAS 
depends on accurate sampling of the gradient artefact 
waveform and timing-alignment of the MR scanner 
and EEG clocks (Mandelkow et al., 2006; Gonsalves et 
al., 2007; Freyer et al., 2009). Benar et al. (2003) also 
indicated that average subtraction processing often 
leaves some remaining artefact inside the signal 
frames. 

The current work presents a novel methodology for 
modelling the gradient artefact waveform to be 
subtracted and restore the EEG signal. Our method 
combines two correction approaches: application of a 
non-linear low-pass filter followed by estimation and 
subtraction of the proposed artefact waveform model 
per se. The non-linear low-pass filter approach, 
modified from Ferreira et al. (2012), makes use of the 
information of high-frequency activity contained in 
the difference between consecutive samples of the 
digital signal. In addition to filter high-frequency 
components, this approach allows higher attenuation 
of low-frequency artefact components as well as less 
distortion in the signal. In turn, the sum of a set of 
sinusoidal waveforms constitutes the basis of the 
proposed artefact waveform model. Our methodology 
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shows to perform effective artefact removal and EEG 
restoration, whose results are very similar to those 
obtained by AAS method. Furthermore, the proposed 
gradient artefact modelling approach shows to predict 
the variability of the artefact waveform over the time, 
which fails to be achieved by the artefact waveform 
averaging. Extensive EPI parameter calculation, data 
segmentation and time alignment are not required as 
well, as shown in the Section Materials and Methods. 

Materials and Methods 

Patients 

Simultaneous EEG and fMRI data have been collected 
for a research focused on epilepsy and post-traumatic 
stress disorder (PTSD) (Van Liempt et al., 2011), jointly 
developed by the department of Psychiatry of 
Universiteit Medisch Centrum Utrecht, the Research 
Centre Military Mental Health Care in the Dutch 
Central Military Hospital in Utrecht, and the 
Department of Research and Development of the 
Kempenhaeghe Epilepsy Center in Heeze, The 
Netherlands. 

Data were recorded from military veterans with PTSD 
which were in mission abroad through the outpatient 
clinic of the Military Mental Health Care. All 
participants were male at the age between 18 and 60 
years. EEG signal recordings from 15 subjects were 
used in order to apply and evaluate the methodology 
we have proposed in this work. 

Protocol for Data Acquisition 

Functional magnetic resonance imaging scanning was 
carried out using a 3 T Scanner (Philips, Eindhoven, 
and The Netherlands) at Kempenhaeghe Epilepsy 
Center. An MRI-compatible 64 channel 
polysomnograph was used to collect one ECG channel, 
two EOG channels, one EMG and 60 EEG channels. 
The EEG electrodes positioning was in accordance 
with the 10 - 20 international system electrodes 
placement. The sampling frequency for EEG signal 
acquisition was 2048 Hz. 

Combined EEG-fMRI data were collected during 45 
minutes, after application of the EEG cap. The subjects 
were scanned using a functional echo-planar imaging 
sequence with 33 transversal slices (thickness 3 mm, 
TE 30 ms, and TR 2500 ms). 

Block Diagram of the Proposed Methodology 
Implementation of our methodology for waveform 
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artefact modelling and EEG restoration was divided 
into four steps, according to the algorithm structure 
block diagram of Fig. 1: 



FIG. 1 PROPOSED ALGORITHM STRUCTURE FOR GRADIENT 
ARTEFACT MODELLING AND EEG RESTORATION 

Each step was developed and applied to the EEG 
recordings in MATLAB environment. 

Peak Identification for ST Estimation 

Estimation of the EPI parameter slice time (ST) was 
performed for determination of the EEG window 
length to be processed as well as the fundamental 
frequencies of the frequency bins (Niazy et al., 2005). 
Observation of the raw data recorded within the MR 
scanning revealed the occurrence of the steep wave 
activity or typical peaks corresponding to the onset of 
each MR slice. Thereby, ST could be estimated by 
evaluation of the time interval between two 
subsequent peaks (Garreffa et al., 2003). 

In order to detect the peaks, initially, a signal 
threshold has been estimated according to the peak 
detection algorithms proposed by Garreffa et al. (2003). 
Considering the data under analysis, the value of ST 
was estimated at 155 ± 1 samples, which corresponded 
to the time interval of 75.68 ± 0.50 ms (Ferreira et al., 
2012 ). 

The EEG window length was set as a multiple of ST. 
This is a required condition during estimation of the 
underlying gradient artefact components within the 
EEG bandwidth, as described below. Hence, an EEG 
window of 32 slices was employed, i.e., 32 x 155 
samples. 

Non-linear Low-pass Filtering Using Signal Slope 
Adaption (SSD) 

To remove high-frequency artefact components, a non- 
linear low-pass filter was developed and utilized 
based upon the difference of consecutive samples of 
the digital signal, the signal slope adaption - SSD 
(Ferreira et al., 2013). Such a filter was implemented by 
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taking into account the idea of estimating the 
variability of the residuals resulting from application 
of the AAS method, as described in Ferreira et al. 
(2012). We would rather use such a filter because, in 
addition to being a zero-phase filter, it provokes less 
distortion in the frames of the filtered signal as well 
causes higher attenuation of the underlying artefact 
components in the EEG bandwidth. These features 
leaded to the minimization of the occurrence of 
artefact residuals in the restored EEG. 

The analysis of the signal slope represented by the 
difference between consecutive samples can be 
employed to detecting high-frequency properties and 
artefacts in digital signals (Scherg, 1982; Barlow, 1983; 
Cluitmans et al., 1993; Van de Velde et al., 1998; 
Ferreira et al., 2012). Cluitmans et al. (1993) and Van 
de Velde et al. (1998) proposed a methodology which 
makes use of such a signal parameter, associated with 
large signal magnitudes, for detection of the high- 
frequency activity related to the muscle artefact in the 
EEG signal. In the same way, observation of the EEG 
recorded within the fMRI scanner reveals that large 
slopes and magnitudes of the signal can be attributed 
to the gradient artefact and its residuals as well 
(Koskinen and Vartiainen, 2009; Ferreira et al., 2012). 

According to Ferreira et al. (2012), the larger slopes 
related to the sharp wave activity of the gradient 
artefact residuals are used to identify whether the EEG 
samples are artefact free or not. In order to perform 
the identification, a slope threshold ( thrs ) is estimated 
in such a way that if the sample has signal slope larger 
than thrs, it is then identified as containing artefact 
interference, thrs can be estimated, for example, taking 
into account the probability distribution of the signal 
slope (Cluitmans et al., 1993; Van de Velde et al., 1998). 
As discussed by Ferreira et al. (2012), it is possible to 
approximate the raw EEG (EEGraw) to the true EEG 
(EEGtme), using the approach proposed by Klados et al. 
(2009). According to this approach, as the value of a 
filter parameter a; is approximated to the optimal 
parameter value a,, the values of the EEGraw would 
tend to the true EEG: 

lim EEG rawi = EEG truei - (1) 

Here, the difference between consecutive samples of 
the signal, diff (EEGraw), has been assumed as the filter 
parameter di, and the slope threshold thrs as the 
optimal value au Taking into account a window whose 


number of points is equal to the length of the raw EEG, 
the maximum absolute value of the difference between 
consecutive samples of the signal EEGraw corresponds 
to the parameter n which is related to diff (EEG raw) by 
the following expression (Ferreira et al., 2013): 

r i = max \diff (EEG raw )| / (2) 

where i is the subscript of the maximum slope within 
diff (EEG raw). 

The two consecutive samples EEGraw,: and EEGraw, :+ i 
associated with n are adapted as follows: 

EEG corcl , = EEG raw , - L ; , 

EEG corcti+1 = EEG mwi+ j + L ; , 

where 

k = P*r i . (4) 

p corresponds to the adaption factor (0 <p < 1) to be 
applied to the parameter ru An appropriate value of p 
was experimentally determined and set to be 0.925. In 
(3), the sign of L: is set positive when EEGraw, :> EEGraw, 
i+i, and vice-versa. The signal EEGraw in equation (3) is 
then replaced by the modified signal EEGcorct, which 
contains the adapted samples EEGcorct, ; and EEGcorct, i+i. 
(2), (3) and (4) are iteratively recalculated until n < thrs. 
The decreasing value of n calculated at each iteration 
ensures the convergence of (3). 

Estimation of thrs has been performed by taking into 
account probability distributions of an artefact free 
EEG interval (EEGtme) and an EEG excerpt containing 
artefacts, both of which picked up from the same EEG 
channel. The true EEG was obtained from the raw 
EEG recorded within the scanner without gradient 
magnetic field variation. It has been observed from the 
probability distributions that more than around 95% of 
the values for the signal slope of the raw EEG, 
diff (EEGraw), lay outside the range 
fV( E EG„„h+ 3 estimated for the signal slope 

of artefact free EEG segments. Further, the probability 
distribution of the slope parameter estimated for the 
true EEG closely resembles Gaussian distributions 
(Cluitmans et al., 1993), in such a way that the 
confidence interval associated with 

Puffc EEG tr „ ! )+3a A y (EE G tru( ; encompasses approximately 

99.5% of the distribution (Papoulis and Pillai, 2002). 
Thereby, estimation of thrs was performed by Eq. (5): 

thrS = /t rfl y fEEGfn] j + 3° diff (EEG <rm > ■ (5) 


42 



Signal Processing Research Volume 2 Issue 2, June 2013 


www. seipub . org/spr 


Raw EEG signal 


1 i i i i i 

466 467 468 469 470 471 472 

Time (s) 

b EEG filtered by the non-linear SSD approach and by the linear low-pass filter 

473 

i i i i i i 

i i I i i i 

Linear filtering 

466 467 468 469 470 471 472 

Time(s) 

C PSD estimate of the signals of (b) 

473 










Linear filtering 



irfiryr W\» 


















Frequency (Hz) 

d PSD of the residual noise from application of SSD approach B e PSD of the residual noise from application of linear filtering 



0 20 40 60 


100 120 140 160 180 200 

Frequency (Hz) 



0 20 40 60 


100 120 140 160 180 200 

Frequency (Hz) 


FIG. 2 COMPARISON BETWEEN THE APPLICATION OF THE SSD NON-LINEAR LOW-PASS FILTER AND A LINEAR LOW-PASS FILTER 


Fig. 2 illustrates how the SSD approach filters the 
high-frequency artefact components. Fig. 2b depicts 
the resulting signals from application of the SSD 
approach and a linear low-pass filter, FIR, 55- 
coefficients, cut-off frequency of 70 FFz, to the raw EEG 
segment of Fig. 2a. This EEG excerpt was selected 
from the recordings of one specific patient, electrode 
position P3. The standard deviation of the subtraction 
between the corrected signals of Fig. 2b is estimated at 
2.1 pV. In Fig. 2c, the PSD estimates of the signals of 
Fig. 2b are depicted. In turn. Figs. 2d and 2e depict the 
residual noise associated with application of both 
filtering approaches. This figure evidences that the 
SSD filter provokes higher attenuation of low- 
frequency artefact components. Frequency 
components corresponding to the harmonics of the 
repetitive slice sequence (multiples of 13.2 Hz) are well 
attenuated in the bandwidth under the filter cut-off 
frequency (Fig. 2d), whereas PSD estimate of Fig. 2e 
shows higher attenuation of such components around 
70 Hz. 

Modelling and Subtraction of Underlying Artefact 
Components 

As underlying artefact components remain in the 
EEGcorct bandwidth after application of the low-pass 
filter described by Eqs. (2) - (5), elimination of such 
components is still necessary in order to restore the 
EEG signal (Allen et al., 2000; Niazy et al., 2005). In 


this respect, the following model was proposed to 
describe the components of the signal EEGcorct: 

K 

EEG corcli = EEG lniei + £ har kJ + n, ( 6 ) 

*=i 

where EEGtme.i is the estimate of the true EEG cleaned 
from gradient artefacts; hark corresponds to the 
fundamental of the k-th artefact harmonic component 
to be removed amongst the total of K harmonic 
components taken into account ( K = 7); and m 
corresponds to the additional noise. Frequency 
components above K = 7 have been neglected since 
they are removed by the SSD low-pass filter. 

The proposed mathematical model for the artefact 
components hark was based upon the analysis of the 
shape and amplitude of average artefact templates 
obtained from application of the AAS method (Allen 
et al., 2000) to the signal EEGcorct. It has been proposed 
that each artefact harmonic component can be 
approximated by a sinusoidal waveform described as: 

har k = A k cos(2nf k t + 0 k ), (7) 

where 

A 2 x Kct,k\ t (8) 

N 

and Scorct.k is the spectral amplitude (FFT) of the k-th 
gradient artefact harmonic component within the 
EEGcorct. N corresponds to the length of the EEGcorct, 
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which must be a multiple of the slice time, ST. This 
requirement allows estimating the value of Ak 
corresponding to the frequency fk which, in turn, is 
provided by Eq. (9): 

fk= k {f' ( 9 ) 

where /s is the EEG sampling rate. Finally, the phase of 
(7) is calculated by. 

6 k =m k y.n, ( 10 ) 

where mk is a number between 0 and 2, estimated in 
such a way that the subtraction of hark from the 
EEGcorci minimizes the value of the spectral amplitude 
of the k-th gradient harmonic artefact component in 
the resulting EEG signal. 

The modelling of the adjacent frequency components 
to the k-th harmonic has been taken into account, 
which occurs in the frequency bins as well. These 
frequency components arise into the region of ±1 Hz 
around the k-th harmonic (Niazy et al., 2005). Thereby, 
the immediate upper and lower frequency 
components inside the bins were inserted into the 
artefact waveform model and subtracted from the 
EEGcorct. 

After subtraction of the underlying artefact 
components from the EEGcorct, a 14-point moving- 
average (MA) filter (Fig. 1) was employed to smooth 


and, thus, additional correction associated with m in (6) 
has been achieved. Hence, the resulting restored EEG 
corresponds to the signal EEGtme. 

Validation 

Our proposed methodology was applied in raw EEG 
excerpts collected from a set of validation EEG 
channels: FP1, F3, F8, T5, P3, Oz, CPI, FC5, AFz, F6, 
C2, TP7, CP4, and POz. The artefact waveform model 
and the EEG restoration achieved by our methodology 
were compared to those obtained by the AAS method 
(Allen et al., 2000). In addition, the reduction of the 
average spectral power associated to the artefact 
frequency bins was assessed, as performed by Niazy et 
al. (2005). 

Results 

Application of our methodology is illustrated in Fig. 3. 
Fig. 3b shows a representative signal EEGcorct resulting 
from the application of the non-linear low-pass filter in 
a raw EEG excerpt which contains strong artefact 
interference (Fig. 3a). This EEG excerpt was selected 
from the recordings of one subject, electrode position 
F8. In order to model and remove the underlying 
artefact components within the EEG bandwidth, Eqs. 
(6) - (10) were then applied to the EEGcorct, resulting in 
the signal depicted in Fig. 3c. The restored EEG signal, 
after application of the MA filter in the signal of Fig. 3c, 
is shown in Fig. 3d. 


Raw EEG signal 



FIG. 3 APPLICATION OF OUR METHODOLOGY TO (A) A RAW EEG SIGNAL WITH STRONG GRADIENT ARTEFACT INTERFERENCE; 
(B) EEGcort, RESULTING FROM APPLICATION OF THE NON-LINEAR LOW-PASS FILTER USING THE SSD APPROACH; (C) EEG 
SIGNAL AFTER SUBTRACTION OF OUR PROPOSED ARTEFACT MODEL; (D) RESTORED EEG SIGNAL 
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a Restored EEG signals C Estimated artefact waveform along the filtered EEG excerpt 



d Estimated average artefact waveforms 





FIG. 4 COMPARISON BETWEEN THE RESULTS OBTAINED BY OUR METHODOLOGY AND THE AAS METHOD: (A) RESTORED EEG 
SIGNALS (B) PSD ESTIMATE OF THE RESTORED SIGNALS; (C) ARTEFACT WAVEFORMS FOR THE ENTIRE EEG EXCERPT; (D) 

AVERAGE ARTEFACT WAVEFORMS (155 SAMPLES) 


Fig. 4 depicts a comparison between the signal 
restorations obtained from the application of our 
methodology and the AAS method in the EEGcorct of 
Fig. 3b. A strong proximity between the restored EEG 
signals can be observed (Figs. 4a and 4b), which was 
confirmed by the median of the standard deviation of 
the difference between the restored signals by both 
methodologies (at 2.3 pV), and the median of the 
estimated averaged artefact waveforms (at 0.6 pV), 
considering the set of EEG validation channels. 
However, unlike an artefact template, the artefact 
waveform modelled by Eqs. (6) - (10) is estimated over 
the EEG excerpt. As noticed in Fig. 4c, such a model 
shows variability over the time. The averaged 
waveform (blue trace) shown in Fig. 4d did result 
from averaging of our proposed model depicted in Fig. 
4c, segmented and interpolated according to the slice- 
timing in the raw EEG signal (as performed by the 
AAS methodology). 

Yet, it has been observed that as the length of the EEG 
excerpt is increased, the standard deviation of the 
difference between the artefact waveforms of Fig. 4c 
tends to decrease and, therefore, such models are 
approximated. According to the artefact correction 
achieved by our methodology, the median value of the 
imaging artefact (Allen et al., 2000) calculated 
considering the validation EEG channels was 2.6 pV 


pk-pk (prior to the application of the MA filter). 

Fig. 5 shows the power spectrum of the artefact 
waveform estimated by our methodology, depicted in 
Fig. 4c, in which the peaks correspond to the spectral 
power attenuated within the artefact frequency bins. 



Frequency (Hz) 

FIG. 5 POWER SPECTRUM OF THE ARTEFACT MODEL 
ESTIMATED BY OUR METHODOLOGY, DEPICTED IN FIG. 4C 

Considering the EEG validation channels, the median 
spectral power attenuation within the frequency bins 
13.2, 26.4, 39.6, 52.9, 66.1, 79.3, and 92.5 Hz achieved 
by the application of the non-linear low-pass filter was 
respectively, 71.4%, 94.5%, 99.8%, 99.8%, 99.9%, 99.9%, 
and 99.9%. Concerning the frequency bands 0-4 Hz, 
4-8 Hz, 8-12 Hz, and 12-24 Hz, the respective 
median spectral power reduction after the subtraction 
of the proposed artefact waveform model was 0.3%, 
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0.6%, 2.3%, and 68%. The higher attenuation in the 
range 12-24 Hz corresponds to the region in which 
the occurrence of removal of components is associated 
with the frequency bins 13.2 and 26.4 Hz. 

Discussion 

Although the artefact resulting from the gradient 
magnetic fields switching has a strong deterministic 
component (Ritter et al., 2010), it is also affected by 
other sources of interference like scanner vibration, 
subject movements during the studies as well as 
inaccurate slice alignment between EEG and fMRI 
clocks. Because of this, the morphology of the artefact 
waveform is not constant over the time (Gonsalves et 
al., 2007; Sun and Hinrichs, 2009). Thus, investigation 
of a more accurate artefact template during the EEG 
restoration by the AAS method (Allen et al., 2000) has 
been often reported in the literature (Mandelkow et al., 
2006; Con calves et al., 2007; Moosmann et al., 2009; 
Koskinen and Vartiainen, 2009). 

The imaging artefact waveform is approximately the 
differential waveform of the MR and RF gradient 
magnetic fields (Brown and Semelka, 2003; Anami et 
al., 2003; Ritter et al., 2010) which are constituted by 
sinusoidal components as well. Thereby, modelling of 
the artefact as sinusoids appears to be a reasonable 
approach. This idea is confirmed by the analysis of the 
artefact template estimated by the AAS method after 
application of the non-linear low-pass filter described 
by Eqs. (2) - (5). 

Hence, the model described by Eqs. (6) - (10) has been 
put forward which performs estimation of the artefact 
waveform by the sum of a set of sinusoidal waveform 
components. According to our approach, the artefacts 
components contained into the frequency bins regions 
(Niazy et al., 2005) are modelled as sinusoids in the 
time-domain (Eq. 7), whereas the respective sinusoidal 
parameters are estimated in the frequency-domain 
(Eqs. (8), (9), and (10)). As observed in Fig. 4e, the 
averaged artefact waveform obtained from our 
methodology is quite similar to that estimated by 
application of the AAS method. In the same way, the 
overall EEG restoration achieved by both 
methodologies is very similar (Fig. 4). Clearly, it 
confirms the hypothesis described by the modelling 
proposal of Eqs. (6) - (10). On the other hand, the 
model calculated by our methodology shows to 
predict the variability of the artefact waveform 
morphology over the time (Fig. 4c), which is not 
observed for the artefact template obtained from the 
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AAS method. The use of an analytical model instead 
an average template could explain this characteristic 
observed within our modelling approach. 

It is also noticed that as the window of the processed 
EEG is increased, the model estimated by our 
methodology approximates to the model estimated by 
the AAS method. This fact suggests that, as our 
proposed model is constituted by the average of 
sinusoids estimated over the EEG interval, it would 
tend to approach to the average model obtained from 
the AAS method for larger EEG excerpts. Furthermore, 
concerning application of the methodology depicted in 
Fig. 1 in EEG excerpts of the other channels, it is 
revealed that the estimated waveform model showed a 
variation from channel to channel, as expected. 

The usage of the non-linear filter described by Eqs. 
(2) - (5) reveals the achievement of a better attenuation 
of the spectral power associated with the low- 
frequency artefact components (Fig. 2) as well as less 
distortion of in the signal frames. These characteristics 
allowed minimizing the occurrence of artefact 
residuals in the restored EEG. Improvements of the 
computational performance of such a filter shall be 
done in future work. Likewise, our methodology 
should be assessed for other types of fMRI scanners. 
As another suggestion for future work, the possibility 
to combine our method with the AAS methodology 
and other artefact correction techniques should be 
investigated as well. 

Conclusions 

In this work, a mathematical model based upon the 
sum of a set of sinusoidal waveforms is proposed for 
gradient artefact correction during continuous EEG- 
fMRI. Estimation and model subtraction for 
restoration of the EEG signal are performed after 
application of a non-linear filter, which makes use of 
the high-frequency information contained in the larger 
artefact signal slopes. 

In addition to the operation of an EEG restoration 
quite similar to that obtained from the AAS method, 
our methodology predicts the artefact waveform 
variability over the time, instead a template with 
constant morphology. 

Such characteristics could enable the usage of our 
method on seeking a more precise artefact waveform. 
Thereby, its combination with AAS itself and other 
correction approaches could improve the quality of the 
restored EEG signal and, thereby, enlarging the range 
of applications of the combined EEG-fMRI technique. 
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